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Abstract 

We  have  measured  the  relative  yields  of  the  cluster  ions,  (NH3)n-iNH2'*’ 
and  (NH3)nH2''',  produced  by  electron  impact  ionization  of  a  supersonic 
molecular  beam  of  neutral  ammonia  clusters,  as  functions  of  cluster  ion  size, 
neutral  beam  stagnation  temperature,  and  electron  impact  energy.  Our 
observation  of  a  magic  number  at  n=7  in  the  distribution  of  (NH3)n-iNH2* 
cluster  ions  is  interpreted  in  terms  of  an  intracluster  bimolecular 
association  reaction  between  a  nascent  NHa""  daughter  ion  and  an  adjacent  NH3 
solvent  molecule,  giving  rise  to  a  protonated  hydrazine  product  ion.  Our 
observation  of  magic  numbers  at  n=5  and  n=8  in  the  (bTi3)nH2‘*'  cluster  ion 
distribution  is  consistent  with  the  production  of  the  NzHe*  Rydberg  radical 
cation  via  an  intracluster  ion-molecule  reaction.  The  dependence  of  the 
yields  of  these  solvated  Rydberg  radical  cations  on  the  stagnation  temperature 
of  the  neutral  molecular  beam  supports  the  notion  that  extensive  solvation 
destabilizes  these  cluster  ions.  Our  experimental  results  are  consistent  with 
recent  theoretical  predictions  concerning  the  structure  and  stability  of  the 
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IntroductiGn 

Nucleation  processes  and  cluster  formation  are  of  fundamental 
significance  in  our  understanding  of  condensed  matter  physics. >-“3  These 
processes  play  a  central  role  in  such  diverse  physical  phenomena  as  crystal 
growth  in  solution  and  at  interfaces,  soot  formation  in  flames,  aerosol 
formation  in  the  upper  atmosphere,  and  particulate  formation  in  the 
interstellar  environment.  Insights  drawn  from  experiment  and  theory  allow  us 
to  answer,  at  least  in  a  qualitative  way,  questions  regarding  the  stability  of 
clusters  as  a  function  of  size  and  structure,  and  the  relative  importance  of 
each  of  the  various  intermolecular  forces  within  these  clusters.  Some  simp]e, 
intuitively  appealing  theoretical  models  have  recently  been  brought  to  bear 
upon  the  problem  of  atomic  cluster  stability  and  its  dependence  on  cluster 
size.  Because  of  the  fruitful  marriage  of  mass  spectrometry  and  various 
techniques  for  the  generation  of  cluster  beams,  it  has  been  possible  in  some 
cases  to  experimentally  verify  theoretical  predictions  of  these  models. 

The  simplest  of  models  is  based  on  geometric  considerations,^  and  is 
applicable  to  clusters  in  which  the  pair  potential  (i.e.,  the  interaction 
between  the  i-th  and  j-th  particles  within  the  cluster)  is  weak  and 
centrosymmetric,  and  the  total  Hamiltonian  is  closely  approximated  by  the  sum 
over  i  and  j  of  the  spatial  integrals  of  the  pair  potential.  Stability  is 
enhcncerj  by  maximizing  contact  between  particles  (via  solid  packing)  while 
minimizing  surface  area;  hence,  the  growing  clusters  attempt  to  assume  a 
^iherical,  or  nearly  spherical,  geometry.  Structures  predicted  to  have 
enhanced  stability  aue  the  solid  Kackay  icosaheora  containing  H  particles: 

Ni  =  (l/3)(10i3  -  15i2  Hi  -  3). 

-  1,  13,  55,  147,  309.  561,...  (1) 
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where  the  index,  i,  can  take  on  any  positive,  non-zero  integer  value.  The 
predominance  of  these  icosahedral  structures  in  the  growth  of  rare  gas  neutral 
clusters  has  been  demonstrated  by  dynamic  computer  simulations  and  gas-phase 
electron  diffraction  spectroscopy,®  as  well  as  cluster  beam  photoionization 
mass  spectroscopy.® 

A  more  sophisticated  model,  the  so-called  .iellium  model. considers 
clusters  not  as  a  collection  of  neutral  particles  held  by  weak  dispersive 
interactions,  but  as  an  ensemble  of  particles  that  can  be  separated  into 
positively  charged  cores  and  valence  electrons.  The  cores  are  assumed  to  be 
smeared  out  into  a  structureless  " jelly"  of  positive  charge,  permeated  by  a 
uniform  cloud  of  itinerant  electrons.  Potential  functions  are  included  in  the 
Hamiltonian  to  account  for  core-core,  core-electron,  and  electron-electron 
interactions.  Since  the  cores  within  a  cluster  are  assumed  to  be  separated  by 
distances  which  are  large  in  comparison  to  their  sizes,  the  core-core 
interactions  can  be  treated  with  a  simple  Coulomb  potential.  The  core- 
electron  potential  must  contain  both  an  attractive  term  (interaction  of 
valence  electrons  with  the  positive  core)  and  a  repulsive  term  (interaction  of 
valence  electrons  with  core  electrons).  The  net  result  is  that  the  core- 
electron  interactions  are  weak.  The  electron-electron  interactions  are 
treated  within  the  context  of  the  local  density  approximation;®  that  is,  the 
interaction  potential  at  r  is  a  function  of  the  local  charge  density  at  r. 

The  potential  must  account  for  electron  exchange  and  correlation.  These 
terms,  as  well  as  the  term  describing  electron  kinetic  energy,  can  be 
described  parametrically  as  functions  of  r«,  the  radius  of  the  spherical 
volume  element  swept  out  by  a  single  gaseous  electron: 


(4/3)Ti(r»Bu>)®  ■  v/N(e) 


(2) 
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where  ao  is  the  Bohr  radius,  v  is  the  volume  of  the  electron  cloud,  and  N(e) 
is  the  total  number  of  electrons  in.  the  cloud.  The  kinetic,  exchange,  and 
correlation  energies  per  electron  can  then  be  expressed  in  terms  of  r».®  The 
total  Hamiltonian  for  the  system  can  then  be  written  as  the  sum  over  all  cores 
of  the  kinetic  and  Coulombic  energies  of  the  cores,  plus  the  sum  over  all 
electrons  of  the  kinetic,  exchange,  and  correlation  energies  of  the  electrons. 
The  energy  eigenstates  of  this  Hamiltonian  are  discrete  and  resemble  those  of 
a  hydrogenic  atom;  however,  the  fact  that  the  resultant  cluster  potential  is 
not  Coulombic  in  form  allows  the  quantum  number  corresponding  to  orbital 
angular  momentum  to  assume  values  greater  than  or  equal  to  the  principal 
quantum  number.  The  spherical  symmetry  of  the  resultant  potential  leads  to 
degeneracies  in  these  eigenstates,  resulting  in  electronic  shell  structure. 
Enhanced  cluster  stability  is  a  consequence  of  the  closing  of  electronic 
shells  and  depends  therefore  on  the  total  number  of  valence  electrons,  rather 
than  upon  core  number  or  structure.  This  model  predicts  shell  closings  for 
N(e)  =  2,  8,  18,  20,  34,  40,  and  so  on.  This  prediction  has  been 
experimentally  verified  by  the  observation  of  "magic  numbers",  or  enhanced 
populations  at  certain  values  of  N(e),  in  the  abundance  spectra  of  neutral  Na 
and  K  clusters, positive  cluster  ions  of  noble  metals, and  negative 
cluster  ions. 5^3 

While  the  approximations  of  the  jellium  model  are  particularly  valid  for 
modeling  metal  clusters,  they  tend  to  break  down  for  the  treatment  of  clusters 
of  molecules  and  non-metal lie  atoms.  Charge  distributions  in  clusters  held 
together  by  covalent  or  hydrogen  bonds  are  not  expected  to  be  uniform,  and  the 
angular  dependence  of  the  core-core  potential  (viz.,  the  constraints  iiposed 
by  orbital  synnetry)  is  expected  to  give  rise  to  preferential  cluster 
geometries.  For  instance,  carbon  clusters  seem  to  grow  as  hollow  polyhedra 


5 


constructed  of  fused  alicyclic  rings.  Enhanced  stability  is  impeirted  when 
ring  strain  is  minimized;  hence,  one  would  expect  polyhedra  constructed  of 
fused  five-  and  six-membered  rings  to  be  especially  stable.  This  contention 
is  experimentally  corroborated  by  the  observation  of  a  magic  number  in  the 
carbon  cluster  ion  distribution  at  N=60,  attributed  to  the  hollow  dodecahedral 
structure,  "buckminsterfullerene" . 5-®  Likewise,  cluster  ions  which  may  form 
intracluster  hydrogen  bonds  are  found  to  display  magic  numbers,  as  evidenced 
by  the  particularly  stable  cluster  ions:  (NH3)4NH4*  (thought  to  have 
tetrahedral  symmetry^®);  and  the  species,  (H20)2oH30*,  (H20)2o(NH3)bH'^  (m=l  to 
6),  and  (H20)27NH4'^  (thought  to  have  ion  clathrate  structures ) .  Although 
fully  elaborated  at  initio  theories  of  the  structure  and  stability  of  large 
molecular  clusters  remain  just  beyond  the  frontiers  of  comDutational 
tractability,  it  is  certainly  clear  that  we  can  gain  qualitative  insights  on 
molecular  cluster  structure  and  stability  from  a  systematic  search  for  magic 
numbers  in  cluster  abundance  spectra.  We  report  in  this  Article  our  results 
for  a  study  of  product  ions  resulting  from  electron  impact  ionization  of  HHa 
clusters,  and  we  rationalize  these  results  in  terms  of  some  intuitive 
structural  models. 

Experimental  Section 

The  experimental  apparatus  has  been  described  in  detail  elsewhere^o  and 
is  illustrated  schematically  in  Figure  1.  Briefly,  the  apparatus  consists  of 
a  differentially  pumped  Campargue  continuous  molecular  beam  source^i  coupled 
to  an  Extrel  C-50  quadrupole  mass  spectrometer.  Anhydrous  NHs  (Linde,  99.99%) 
is  introduced  into  the  stagnation  volume  of  the  apparatus  in  order  to  maintain 
a  stagnation  pressure  of  1  atm.  Neutral  clusters  fora  in  the  free- jet 
expansion  created  as  the  high-pressure  gas  expands  adiabatically  through  a  250 
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pm  diameter  nozzle  into  the  source  chamber.  The  neutral  cluster  population 
distribution  can  be  regulated  to  a  great  extent  by  changing  the  nozzle 
temperature.  In  these  experiments,  the  nozzle  was  typically  maintained 
(within  ±0.1  K)  at  temperatures  within  the  range  of  253  to  313  R  by  a 
recirculating  chiller.  Pressure  in  the  source  chamber  is  maintained  at  about 
50  mTorr  during  operation  of  the  cluster  beam  by  a  1000  m^/hr  roots  blower. 

The  cluster  beam  passes  from  the  source  chamber,  through  a  0.5  mm  diameter 
skimmer  (about  30  nozzle  diameters  downstream  from  the  start  of  the 
expansion),  and  into  a  colli mat ion  chamber.  This  chamber  is  pumped  by  a  360 
1/sec  turbopump  and  is  maintained  at  10“®  Torr  during  beam  operation.  The 
beam  passes  out  of  the  collimation  chamber  through  a  second  0.5  mm  diameter 
skimmer  and  into  the  mass  spectrometer  chamber.  This  chamber  is  also  pumped 
by  a  360  1/sec  turbopump  and  is  maintained  at  5  x  lO-"^  Torr  during  beam 
operation. 

The  cluster  beam  enters  the  ion  source  and  quadrupole  mass  filter  in  an 
axial  configuration.  The  mass  filter  is  typically  operated  at  better  than  0.3 
amu  mass  resolution  up  to  about  1200  ainu,22  although  we  can  select  the  best 
compromise  between  resolution  and  sensitivity,  depending  on  the  demands  of  a 
particular  experiment.  Emission  current  in  the  source  is  maintained  at  1  mA 
by  current  regulation  at  fixed  electron  energy.  In  these  experiments,  mass 
spectra  were  collected  at  various  electron  kinetic  enersfies  ranging  from  about 
40  to  100  eV.  Cluster  ions  and  fragments  formed  in  the  ion  source  are 
filtered  by  the  quadrupole  and  imaged  onto  an  off-axis  ehanneltrnn.  The 
amplified  signal  from  the  channeltron  is  then  averaged  by  a  digital  storage 
oscilloscope. 

In  order  to  assure  the  validity  of  the  mass  spectral  data,  we  regularly 
calibrated  the  mass  scale  and  the  sensitivity  of  the  quadrupole  mass 
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spectrometer.  In  the  range  below  500  amu,  we  calibrated  the  instrument 
against  the  mass  spectrum  of  perf luorotributylamine;  beyond  500  amu,  we  used 
an  argon  cluster  beam  as  a  calibration  standard. 

Results  and  Discussicn 

Ion  Intensities.  The  high-resolution  electron  impact  mass  spectrum  of 
NHa  clusters  can  be  described  as  an  exponentially  decaying  progression  of 
quartets.  In  the  limit  of  low  electron  kinetic  energy,  this  exponential  decay 
is  primarily  a  reflection  of  the  population  distribution  in  the  neutral  beam. 

A  portion  of  the  mass  spectrum,  which  is  illustrated  in  Figure  2,  shows  the 
quartet  corresponding  to  cluster  ions  containing  25  nitrogen  atoms.  The  peak 
labeled  P  represents  the  unreacted  parent  peEJt,23  while  the  peak  labeled  1 
represents  the  solvated  product  of  the  intracluster  ion-molecule  reaction 
corresponding  to  the  bimolecular  analog: 24 

NHa*  +  NHa  — >  NH4-^  +  NHa  ,  A//*  =  -1.1  eV  (Z) 

Peak  a  corresponds  to  cluster  inns  with  the  emoirinal  formula.  fNH5»')24NH2'^ . 

The  nascent  ions  can  be  thought  of  as  solvated  daughter  ions,  resulting  from 
unimolecular  fragmentation  of  the  solvated  NHa*  parent.  One  might  suspect 
this  daughter  ion  to  undergo  subsequent  intracluster  ion-molecule  reactions  on 
the  timescale  of  the  experiment  (vide  infra).  Peak  a'  corresponds  to  a 
cluster  ion  with  the  empirical  formula,  (NHa)2sH2*.  It  has  been  suggested 
that  this  cluster  ion  represents  the  solvated  product  of  an  unusual 
intracluster  ion-molecule  reaction  having  no  known  gas-phase  bimolecular 
analog. 2C  Although  a  small  part  of  the  intensitv  of  this  peak  imv  he 
attributed  to  the  presence  of  H2O  in  the  clusters2B  (due  to  water 
contamination  in  the  stagnation  chamber),  results  from  t.he  experiments  herein, 
as  well  as  results  from  previous  experiments  in  our  laboratory  on  doubly 
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charged  NHa  clusters, 2°  support  the  identity  of  this  cluster  ion  as  an  unusual 
nitrogen-containing  species  formed  via  an  intracluster  ion-molecule  reaction. 

A  mechanism  accounting  for  the  mass  spectral  features  described  above  is 
shown  schematically  in  Figure  3.  We  now  wish  to  consider  the  relative 
abundance  of  these  cluster  ions  as  a  function  of  cluster  size.  Castleman  and 
co-workersi®  have  examined  the  population  distribution  of  the  cluster  ions, 
(NH3)n-iNH4'^,  and  have  observed  a  “magic  number"  at  n=5,  corresponding  to  the 
particularly  stable  tetrahedral  cluster  ion,  (NH3)4NH4"^.  We  discuss  below  the 
appearance  of  "magic  numbers"  in  the  population  distributions  of  the  ions, 
and  (NHslnHz'*’,  and  interpret  these  findings  in  terms  of 
structural  models. 

The  (NH3)n-iMH2r^  Cluster  Ions.  A  plot  of  the  relative  intensities  of 
the  (NHaln-iNHz'*’  ions  as  a  function  of  cluster  size  n  is  shown  in  Figure  4.  A 
magic  number,  corresponding  to  enhancement,  in  the  inn  Riimal .  is  nhserved  for 
the  cluster  ion  of  size  n=7.  This  enhancement  at  n=7,  relative  to  ion  yield 
at  n  >  11,  is  independent  of  electron  impact  energy,  although  absolute  ion 
yields  are  observed  to  increase  as  a  function  of  electron  impact  energy.  As 
shown  in  Figure  5,  this  relative  enhancement  appears  to  increase  somewhat  as 
the  nozzle  temperature  is  increased.  This  is  presumably  due  to  a  more  rapid 
exponential  decay  of  the  neutral  cluster  beam  population  distribution  (i.e.,  a 
smaller  number  of  larger  clusters)  under  conditions  of  higher  stagnation 
temperatures . 

We  can  account  for  the  observed  size  dependence  of  the  (NHaln-iNHz* 
cluster  ion  yield  if  we  assume  that  an  associative  ion-molecule  reaction 
occurs  between  the  nascent  NHz'*’  ion  and  an  adjacent  NHa  solvent  molecule 
within  the  cluster: 


NHg*  +  NHa - >  NgHs- 


(4) 


From  available  thermochemical  data.^v.ze  estimate  that  this  reaction  is 
exothermic  by  4.52  eV.  A  similar  (albeit  endothermic)  bimolecular  reaction  is 
known  to  occur  between  NHa*  and  NHs: 

NHa*  +  HHa - >  N2H5*  +  H  ,  tff"  -  1.12  eV  (5) 

at  a  rate  which  is  about  0.1%  cf  gas-kinetic. 2® *3°  We  expect  the  barrier  for 
reaction  (4)  to  be  smaller  than  that  for  reaction  (5)  for  the  following 
reasons.  First,  within  the  context  of  the  Hammond  postulate, 3^  the  barrier  to 
reaction  (4)  is  expected  to  occur  earlier  along  the  reaction  coordinate,  the 
potential  surface  is  expected  to  be  more  attractive,  and  the  transition  state 
is  expected  to  more  closely  resemble  the  reactants.  Second,  the  activation 
energy  for  reaction  (4)  is  expected  to  be  lower  than  that  of  reaction  (5), 
since  the  transition  state  for  (5)  must  contain  at  least  one  highly  extended 
N-H  bond  (correlating  with  eventual  loss  of  the  product  H  atom).  By  analogy 
with  other  associative  bimolecular  reactions,  we  would  expect  the  activation 
energy  for  reaction  (4)  to  be  small,  if  not  negligible,  and  determined  to  a 
great  extent  by  the  dipole-dipole  interaction . 32  Third,  the  symmetry 
contribution  to  the  barrier  for  (5)  is  probably  larger  than  that  for  (4),  due 
to  the  N-H  bond  distortion  (and  subsequent  reduction  in  symmetry  number)  in 
the  transition  state  for  (5).  We  might  therefore  expect,  based  on  differences 
in  activation  energies  and  Arrhienius  A  factors,  that  the  associative  reaction 
(4)  will  proceed  at  a  rate  significantly  greater  than  0.1%  of  the  gas-kinetic 
rate. 

While  the  "naked"  (unclustered)  product  of  the  highly  exothermic 
associative  reaction  (4)  would  have  sufficient  internal  energy  to  undergo 
subsequent  H-H  bond  cleavage,  it  would  certainly  be  stabilized  through 
solvation  by  additional  NHs  molecules.  In  fact,  if  we  hydrogen-bond  five  NHs 
solvent  molecules  to  the  five  H  atoms  of  the  oroduct  ion.  we  end  up  with 
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an  ion  having  a  completed  solvation  shell  with  the  empirical  formula, 
(NH3)8NH2"^,  as  illustrated  in  Figure  6.  This  cluster  ion  has  a  total  of  seven 
N  atoms  and  accourits  for  our  observation  of  a  magic  number,  n=7. 

The  (HH3)nHz^  Cluster  lens.  A  olot  of  the  intenKitiRs  nf  the  fNH3)nH2‘^ 
ions  as  a  function  of  cluster  size  n  is  shown  in  Figure  7.  Two  magic  numbers 
are  observed  in  the  cluster  ion  population  distribution,  corresoonding  to 
enhanced  stability  of  the  n=5  and  n=8  cluster  ions.  The  relative  enhancements 
at  n=5  and  n=8  (compared  to  cluster  ions  of  size  n  ^  11)  are  independent  of 
electron  impact  energy,  although  absolute  ion  yields  are  again  observed  to 
increase  at  higher  energy.  These  enhancement.^  noDenr  tn  increase  dramatically 
as  the  stagnation  temperature  is  increased  from  253  to  273  K,  as  shown  in 
Figure  8.  The  increase  in  relative  enhancement  as  a  function  nf  increasing 
stagnation  temperature  is  more  gradual  at  temperatures  above  273  K.  This 
temperature  dependence  will  be  discussed  in  greater  detail  helnw. 

One  might  argue  that  we  are  not  really  observing  cluster  ions  of  the 
form  (NH3)nH2*,  but  rather  the  isobaric  roecies  .zb  Such 

heteroclusters  can  be  thought  of  as  solvated  NH4'*'  ions  in  which  one  of  the  NHa 
solvent  molecules  has  been  replaced  by  an  H2O  molecule.  These  heteroclusters 
would  be  expected  to  have  structures  analogous  to  those  of  the  (NH3)n-iNH4'^ 
cluster  ions,  and  would  therefore  display  magic  numbers  in  the  population 
distriPution  at  n=5  (packed  tetrahedron )i®  and  n=9  (packed  dodecahedron ), 33  as 
observed  for  solvated  NH4*  cluster  ions.  However,  we  do  not  observe  magic 
numbers  at  these  n  values,  but  rather  at  n=5  and  n=8.  Therefore,  we  believe 
that  what  we  are  observing  is  the  series  of  cluster  ions,  (NH3>nH2'^,  and  not 
isobaric  interferences  from  HzO-containing  heteroc lusters. 

Such  a  molecular  formula  implies  some  interesting  structural 
possibilities.  Garvey  and  Bernstein  originally  suggestedzs  that  these  cluster 
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ions  night  arise  from  intracluster  ion-molecule  reactions  which  give  rise  to  a 
solvated  a  hypervalent  species  with  trigonal  bipyramidal  symmetry.  Such 

a  structure  is  consistent  with  the  results  of  ab  initio  calculations,®* 
although  a  minimum  in  the  hypersurface  is  achieved  only  if  the  system  is 
"Rydbergized"  .®®  We  might  expect  that  an  NHs"*  ion  solvated  by  five  NHa 
solvent  molecules  would  show  enhanced  stability,  by  analogy  with  (NH3).4NH4*, 
and  that  the  population  distribution  of  (NH3)n-iNH6*  cluster  ions  would 
therefore  display  a  magic  number  at  n=6.  However,  we  do  not  observe  a  magic 
number  at  n=6,  and  we  cannot  account  for  our  observed  magic  numbers  on  the 
basis  of  solvation  of  a  trigonal  bipyramidal  species.  Recent  ab  initio 
calculations  at  the  unrestricted  Hartree-Fock  level®®  suggest  that  the 
unsolvated  ion  is  actually  a  Rydberg  radical  ion,  NaHe*,  whose  predicted 
structure  is  shown  in  Figure  9.  One  might  expect  that  an  NzHe*  ion  solvated 
(via  hydrogen  bonding)  at  each  of  its  six  available  H  sites  would  display 
enhanced  stability  with  respect  to  species  of  different  degrees  of  solvation. 
The  population  distribution  of  (NHaln-zCHzHe)*  cluster  ions  would  therefore 
show  a  magic  number  at  n=8;  this  is  precisely  what  we  observe  experimentally. 
Ojr  observation  of  magic  numbers  at  n=5  and  n=8  can  be  accounted  for  on  the 
basis  of  the  two  structures  shown  in  Figure  10.  Apparently,  enhanced 
stability  is  imparted  even  when  only  one  of  the  ion's  nitrogen  centers  is 
coopletely  solvated;  that  is,  when  the  solvation  shell  surrounding  the  NaHe'^ 
ion  is  only  half-full. 

Kassab  et  al.  suggest,®®  by  analogy  with  the  NH4  Rydberg  radical,®*^  that 
extensive  solvation  of  HzHa"'  may  actually  destabilize  the  radical  cation 
cluster.  Hence,  one  might  expect  to  observe  (NH3)n-2N2Her*'  cluster  ions  only 
for  small  to  moderate  values  of  n.  If  one  were  to  prepare  these  cluster  ions 
by  electron  impact  ionization  of  a  molecular  beam  of  neutral  clusters  of  NHs 
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(as  we  have  done),  one  might  then  expect  to  observe  greater  yields  of 
(NH3)n-2N2H0"*-  when  the  neutral  beam  has  an  enriched  population  of  smaller 
clusters.  One  way  to  attain  such  an  enrichment  is  to  operate  the  molecular 
beam  expansion  at  elevated  stagnation  temperatures.  In  fact,  we  do  observe  an 
increase  in  the  yields  of  these  cluster  ions  as  the  stagnation  temperature  is 
increased,  as  shown  in  Figure  8.  This  observation  represents  experimental 
confirmation  of  the  theoretical  predictions  of  Kassab  et  al.^^ 

Conclusions 

He  have  reported  previously  unobserved  magic  numbers  in  the  population 
distributions  of  the  cluster  ions,  (NH3)n-iNH2'*-  and  (NH3)nH2*,  produced  via 
electron  impact  ionization  of  a  supersonic  molecular  beam  of  neutral  (NH3)n+i 

i 

clusters.  We  take  the  observation  of  a  magic  number  of  n=7  in  the  solvated 
NH2*  cluster  ion  distribution  as  evidence  of  an  associative  ion-molecule 
reaction  between  an  daughter  ion  and  an  NHa  solvent  molecule,  giving  rise 
to  a  protonated  hydrazine  molecule  within  an  NHa  solvent  shell.  The  magic 
numbers  we  observe  in  the  (NHalnHz'*’  cluster  ion  population  distribution  can  be 
explained  in  terms  of  a  solvated  NzHe*  Rydberg  radical  cation,  and  is 
inconsistent  with  production  of  solvated  heterocluster  ions  of  the  form, 
(HzOXNHsln-^H^*.  Our  results  are  ccxisistent  with  the  ab  initio  calculations 
of  Kassab  et  al.  on  the  structure  and  stability  of  the  NzHe'*'  Rydberg  radical 
cation.  It  is  hoped  that  continuing  experimental  efforts,  in  our  own 
laboratory  as  well  as  those  of  our  colleagues,  will  stimulate  the  further 
elaboration  and  refinement  of  theoretical  models  dealing  with  the  electronic 
structure  of  clusters  and  cluster  ions. 


13 


Acknovled^ent  is  made  to  the  Office  of  Naval  Research,  and  to  the 
donors  of  the  Petroleum  Research  Fund,  administered  by  the  American  Chemical 
Society,  for  generous  support  of  this  research. 


REFERENCES 


(1)  Hoare,  M.R.  Adv.  Chem.  Phys.  1979,  40,  49. 

(2)  Mftrk,  T.D.;  Castleman.  A.W..  Jr.  Adv.  At.  Hnl .  P/tvr.  19B5.  20,  65. 

(3)  Mfixk,  T.D.  Jnt.  J.  Mass  Soectrom.  Inn  Prnn.  lfiR7.  7f).  1 . 

(4)  Mackay.  A.L.  Anta  Crvst.  19R2.  7.*).  91 R. 

(5)  Farges,  J.;  deFeraudv.  M.F.:  Rnniilf..  R.  :  Tnrph#>t.  0.  Tntamational 
Meeting  on  Small  Particles  and  Inorganic  Clusters,  Lyon-Villeurbanne, 
September  1976.  Journal  de  Physique  Colloaue  C-2.  1977.  d.  47. 

(6)  Schriver,  K.E. :  Hahn.  M.Y.;  Persson,  J.L.;  LaVilla,  M.E.;  Whetten, 
R.L.  J.  Phys.  Chem.  1989,  93,  2869. 

(7)  Cohen,  H.L.;  Chou,  M.Y.;  Knight.  W.D.;  de  Heer,  W.A.  J.  Phys.  Chem. 
1967,  91,  3141. 

(8)  Kohn,  W.;  Sham,  L.J.  Phys.  Rev.  1965,  140,  A1333. 

(9)  Kittel,  C.  Introduction  to  Solid  State  Physics,  6th  ed.;  Wiley:  Hew 
York,  1986. 

(10)  Knight,  W.D.;  Clemenger,  R.;  de  Heer,  W.A.;  Saunders,  W.A.  Solid 
State  Common.  1885,  53,  445. 

(11)  Katakuse,  I.;  Ichihara,  I.;  Fujita,  Y. ;  Matsuo,  T.;  Sakurai,  T. ; 
Matsuda,  H.  Int,  J.  Mass  Spectrom.  Ion  Proc.  1965,  67,  229. 

(12)  Begemann,  W.;  Meiwes-Broer,  K.H.;  Lutz,  H.O.  Z.  Phys.  D  1966,  3, 

109. 

(13)  Horti^,  G.;  Mdller,  M.  Z.  Phys.  1969,  221,  119. 

(14)  Joyes,  P.;  Sudraud,  P.  Surf.  Sci.  1985,  156,  451. 

(15)  Kroto,  H.W.;  Heath,  J.R.;  O'Brien,  S.C.;  Curl,  R.F.;  Smalley,  R.E. 
Nature  1965,  318,  162. 

(16)  Echt,  0.;  Morgan,  S.;  Dao,  P.D.;  Stanley,  R.J.;  Castleman,  A.W., 
Jr.  Ber.  Bunsenges.  Phys.  Chem.  1864,  88,  217. 


(17)  Kassner,  J.L.;  Hagen,  D.E.  J.  Chem.  Phys.  1976,  64,  1860. 

(18)  Shinohara,  H.;  Nagashijna,  U.;  Nishi,  N.  Chem.  Phys.  Lett.  1964, 


15 


111,  511. 

(19)  Shinohara,  H.;  Nagashima,  U.;  Tanaka,  H.;  Nishi,  N.  J.  Chem.  Phys. 
1965,  83,  4183. 

(20)  Coolbaugh,  M.T.;  Peifer,  H.R.;  Garvey,  J.F.  Chem.  Phys.  Lett.  1989, 
156,  19. 

(21)  Campargue,  R.  J.  Phys.  Chem.  1984,  88,  4466. 

(22)  Peifer,  W.R.;  Coolbaugh,  M.T.;  Gsurvey,  J.F.  J.  Phys.  Chem.  1969, 

93,  in  press . 

(23)  It  is  also  possible  that  this  "parent"  peak  represents,  at  least  in 
part,  the  solvated  product  ion  of  Reaction  (3);  that  is,  (NH3)n-2NH2NH4'*' .  For 
larger  clusters,  reaction  exothennicity  nay  be  carried  away  either  by  loss  of 
the  neutral  NH2  product,  or  by  evaporation  of  a  surface  NHs  solvent  molecule. 
In  the  latter  case,  it  is  not  possible  to  distinguish  the  unreacted  parent 
cluster  ion  from  the  cluster  ion  which  has  undergone  a  bimolecular  reaction. 
See,  for  example:  Morgan,  S.;  Castleman,  A.W.,  Jr.  J.  Am.  Chem.  Soc.  1967, 

109,  2867. 

(24)  Hogg,  A.M.;  Haynes,  R.H.;  Rebarle,  P.  J.  Am.  Chem.  Soc.  1966,  86, 

28. 

(25)  Gsirvey,  J.F.;  Bernstein,  R.B.  Chem.  Phys.  Lett.  1968,  143,  13. 

(26)  Shinohara,  H.;  Nishi,  N.  Chem.  Phys.  Lett.  1967,  141,  292. 

(27)  Rosenstock,  H.M.;  Draxl,  K. ;  Steiner,  B.W.;  Herron,  J.T.  J.  Phys. 
Chem.  Ref.  Data  1977,  6,  Suppl.  Monograph  1. 

(28)  Meot-Ner,  M.;  Nelsen,  S.F.;  Willi,  M.R.;  Frigo,  T.B.  J.  Am.  Chem. 


Soc.  1964,  106,  7384. 


16 


(29)  Deivish,  6.A.W.;  A.;  Giardini-CJuidoni,  A.;  Volpi,  G.G.  J. 

Chea.  Phys.  1963,  39,  1599. 

(30)  Melton,  C.E.  J.  Chea.  Phys.  1966,  45,  4414. 

(31)  Hammond,  G.S.  J.  Am.  Chea.  Soc.  1955,  77,  334. 

(32)  Benson,  S.W.  Theraocheaical  Kinetics,  2nd  ed.;  Wiley:  New  York, 

1976. 

(33)  In  addition  to  the  magic  number,  n=5,  observed  by  Echt  et  al.  for 
the  (NH3)n-iNH4'^  clusters  ions,  we  observe  an  additional  magic  number  at  n=9. 
We  attribute  this  observation  to  the  packing  of  four  additional  NHa  molecules 
along  the  four  faces  of  the  tetrahedron  of  (NH3)4NH4'^,  although  we  cannot  rule 
out  the  possibility  of  an  alternative  packing  configuration. 

(34)  Bugaets,  O.P.;  Zhogolev,  D.A.  Chea.  Phys.  Lett.  1977,  45,  462. 

(35)  Mulliken,  R.S.  Accts.  Chea.  Res.  1976,  9,  7. 

(36)  Kassab,  E.;  Fouquet,  J.;  Evleth,  E.M.  Chea.  Phys.  Lett.  1988,  153, 

522. 

(33)  Kassab,  E.;  Evleth,  E.M.  J.  Aa.  Chea.  Soc.  1987,  109,  1653. 


Schematic  representation  of  the  differentially 


pumped  cluster  beam  apparatus  and  quadrupole  mass 
spectrometer.  The  temperature  of  the  nozzle  in  the 
stagnation  region  is  regulated  by  a  circulating  chiller. 
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A  representative  portion  of  the  70  eV  electron 


impact  mass  spectrum  of  the  ammonia  cluster  beam.  This 
portion  of  the  spectrum  corresponds  to  the  four  different 
cluster  ion  species  containing  25  nitrogen  atoms  each.  The 
abscissa  is  scaled  according  to  the  nominal  masses  for  H  and 
N  of  1  and  14  amu ,  respectively. 
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Schematic  representation  of  the  dynamical 


processes  leading  to  formation  of  the  various  solvated 
product  ions  from  the  solvated  parent.  Differences  in  the 
shading  of  the  indicated  peaks  is  intended  to  indicate  that 
the  four  peaks  do  not  belong  to  the  same  family  of  n-mers. 
Peaks  P  and  a  are  members  of  the  family  of  (n+l)-mers,  whil 
peaks  1  and  a'  are  members  of  the  family  of  n-mers. 
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A  plot  of  the  intensities  of  the  (NH3)n-iNH2 


cluster  ions  versus  cluster  size,  n.  The  stagnation 
temperature  was  held  at  273  K  during  the  collection  of  all 
mass  spectra  represented  here.  Note  the  magic  number  at 
n  =  7. 
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A  plot  of  the  intensities  of  the  (NHo )n-iNH2‘^ 


cluster  ions  versus  stagnation  temperature.  The  mass 
spectra  represented  here  were  all  collected  at  an  electron 
impact  energy  of  70  eV. 
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Figure  6.  Proposed  structure  for  the  (NH3)8NH2'^  cluster 
ion.  This  species  is  the  nost  prevalent  of  all  cluster  i 
in  the  series,  (NHa )n-iNH2‘^ ,  and  is  believed  to  be  a 
protonated  hydrazine  molecule  surrounded  by  one  complete 
solvation  shell. 
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Figure  7.  A  plot  of  the  intensities  of  the  (NH3)nH2‘^ 
cluster  ions  versus  cluster  size.  n.  The  stagnation 
temperature  was  held  at  273  R  during  the  collection  of  all 
mass  spectra  represented  here.  Mote  the  magic  numbers  at 


n=5  and  n=8. 
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Figure  8.  A  plot  of  the  intensities  of  the  (NH3)nH2* 
cluster  ions  versus  stagnation  temperature.  The  mass 


spectra  represented  here  were  all  collected  at  an  electron 
impact  energy  of  70  eV.  Note  that  the  enhancement  of  the 
"magic"  ion  signals,  relative  to  signals  from  the  species 
larger  than  n=ll,  appears  to  increase  dramatically  as  the 

stagnation  temperature  is  increased  from  253  K  to  273  K. 

The  degree  of  enhancement  increases  more  gradually  as  the 

stagnation  temperature  is  increased  from  273  R  to  313  K. 
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Figure  9. 
et  al.^^ 


Proposed  structure  of  the  N2Hs'^  ion,  after  Kassab 
This  species  is  predicted  to  be  a  Rydberg  radical 


cation . 


Figure  10.  Proposed  structures  for:  a)  the  (NH3)3N2Hb'^ 
cluster  ion,  and  b)  the  (NHaleNaHe'^  cluster  ion.  These 
species  represent  the  nost  stable  solvated  forms  of  the 
Rydberg  radical  cation,  Apparently,  complete 

solvation  of  one  end  of  this  radical  ion  is  sufficient  to 
impart  enhanced  stability. 


